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Abstract 

We solve the escape problem for the Heston random diffusion model. We obtain exact expressions 

for the survival probability (which ammounts to solving the complete escape problem) as well as 

66' 

^^ ■ for the mean exit time. We also average the volatility in order to work out the problem for the 

o: 

CN . return alone regardless volatility. We look over these results in terms of the dimensionless normal 

level of volatility -a ratio of the three parameters that appear in the Heston model- and analyze 

C^^- ' their form in several assymptotic limits. Thus, for instance, we show that the mean exit time grows 

quadratically with large spans while for small spans the growth is systematically slower depending 
on the value of the normal level. We compare our results with those of the Wiener process and 
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l^\ [ show that the assumption of stochastic volatility, in an apparent paradoxical way, increases survival 
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and prolongs the escape time. 
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I. INTRODUCTION 

Models of financial dynamics based on two-dimensional diffusion processes, known as 
stochastic volatility (SV) models [l[, are being widely accepted as a reasonable explanation 
for many empirical observations collected under the name of "stylized facts" J2]. In such 
models the volatility, that is, the standard deviation of returns, originally thought to be 
a constant, is a random process coupled with the return so that they both form a two- 
dimensional diffusion process governed by a pair of Langevin equations l|. 

Volatility is nowadays a key magnitude in any financial setting. It is the backbone of 
many financial products that are designed to cover investors' risk. Extreme values associated 
with volatility have thus a special meaning, as they do in physics and natural sciences 
where escape problems in noisy environments such as Kramers problem are of the utmost 



importance 



\i. 



In a recent paper we have addressed a partial aspect of the problem: that of extreme times 
for the volatility regardless the value of the return [a] . Now we want to address the overall 
escape problem associated with both return and volatility. This is certainly a more difficult 
task because the return strongly depends on volatility while, in the standard approach to 
SV models, the latter is supposed to be independent of the former. 

We are thus left with a two-dimensional escape problem which is always quite involved. 
The situation is similar to that of the unbounded Brownian particle where the extreme- value 

' the particle is relatively easy to handle while that of its position 
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problem for the velocity_p: 
is much more intricate 

The extreme-time problem of the return has been addressed, to our knowledge, in only 
few works. We refer the reader to our recent work on the subject (9|, [lOj, although it is based 
on the continuous-time random walk technique which is an entirely different frame, yet with 
a different scope, than that of SV models. Within the setting of the latter we are only aware 
of the recent work of Bonano et al 
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12| in which an approach to the escape-time problem 



is addressed through simulations of the Heston model. 

In this paper we study the complete extreme value problem of one particular SV model: 
the Heston model [ij] . Different S V models basically differ in the way the volatility depends 
on the underlying noise governing its dynamics. The Heston model has the benefit, over 
other SV models, of allowing exact analytical developments. This is the case of its unre- 



stricted {i.e., barrier-free) probability density function which was obtained by Yakovenko 
and Dragulescu few years ago [l3|. Herein we will obtain not only the exact expression 
of the the mean escape time (MET) but the exact survival probability as well. Being the 
knowledge of the latter equivalent to solving the entire escape problem. 

The paper is organized as follows. In Sec. [Ill we present the Heston model and obtain 
the complete solution to the escape problem. In Sec. IIIII we evaluate the mean escape time 
and analyze its behavior for high and low volatility. In Sec. [IV]we average out the volatility 
assuming it has reached the stationary state. This allows us to get exact expressions for the 
survival probability and the mean escape time of the return alone. Conclusions are drawn 
in Sec. |V] and some more technical details are in appendices. 

II. THE HESTON MODEL AND THE SURVIVAL PROBABILITY 

Let P{t) be a speculative price or the value of a financial index. We define the zero-mean 
return X{t) through the stochastic differential: 



where (■) denotes the average. In terms of X{t) the Heston model l4] is a two-dimensional 
diffusion process {X(t), Y{t)) described by the following pair of stochastic differential equa- 
tions 



dX{t) = ,/Yit)dW,{t), (2) 

dY{t) = -a [Y{t) - m^] dt + k^/Y{{}dW2{t), (3) 

where Wi(t) are Wiener processes, i.e. dWi(t) = C,i{t)dt {i = 1,2), where C,i{t) are zero-mean 
Gaussian white noises with {^i(t)^i(t')) = 6ij6(t — t') [15|]- Note that in this particular model 
the volatility is 



c^(t) = V^, (4) 

i.e., Y{t) is the variance of return although, as long as no confusion arises, we will use the 
term "volatility variable" or just "volatility" for the random process Y{t). In Eq. ([3]) the 
parameter m is the so-called normal level of volatility, a > is related to the "reverting 
force" toward the normal level (see below) and k, sometimes referred to as the "vol-of-vol" , 
measures the fluctuations of the volatility. 



In the context of biological diffusion problems the process Y(t) described by Eq. (^ was 



proposed many years ago by Feller 16j who, among other properties, proved that Y{t) is 
always positive so that the volatility, Eq. (jlj), is real, positive and well defined. This feature 
along with a non-negligible (and exponential) autocorrelation with characteristic time 1/a 
makes the process very appealing from the perspective of mathematical finance. 

In 1985 Cox, Ingersoll and Ross [iTJ introduced the same dynamics in connection with 
interest rates of bonds. Almost a decade later and aiming to provide a more realistic price 
for options, Heston [ij] undertook the same dynamics but for the diffusion coefficient of 
financial price fluctuations as is precisely shown in Eqs. (El)- ([3]). 

The resulting process has become quite popular among financial practitioners who want 
to include the effect of volatility changes in option pricing. Part of this success is due to 
the easy interpretation of the parameters. As mentioned, 1/a provides the typical time 
that the volatility needs to reach the stationary state (the stationary density is the Gamma 
distribution, see Sec. [IV|) . For this reason, a can also be interpreted as the strength of 
the reverting force that ties the process Y{t) to its normal level m^, the latter being the 
mean value of Y{t) in the stationary state. The magnitude of the volatility fluctuations is 
provided by k which like a and m^ have all units of l/(time). 

Our main interest is the escape problem associated with the Heston model. To this 
end, let us denote by S{x,y,t) the probability that the zero-mean return X{t), starting 
at X{0) = X with volatility Y{0) = y, is at time t inside the interval (— L/2,L/2) without 
having ever left it during previous times. In other words, S{x, y, t) is the survival probability 
(SP) for the joint process {X(t)^ Y{t)) to be at time t inside the strip 

-L/2 < X{t) < L/2, < Y{t) < cx), 

with X(0) = X and F(0) = y. 

The SP obeys the following backward Fokker-Planck equation [1 

dS , 2^^'^ 1,2 Q'^S 1 d'^S 

m = -"^^ - "^ ) a^ + 2^ ^9^ + 2^9^' (^) 

with initial and boundary conditions respectively given by 

S{x,y,0) = l, S{±L/2,y,t) = 0. (6) 

This problem can be solved by means of Fourier series. Indeed the boundary conditions. 



S{±L/2, y, t) = 0, lead us to look for a solution of the form 

oo 

S{x, y, t) = ^ Sn{y, t) COS [{2n + 1) TTx/L] , (7) 



Sn{y,t) = — / S'(x, y, t) cos [(2ri + 1) vrx/L] (ix. (8) 

-^ i-L/2 



n=0 

where the Fourier coefficients Sn{y,t) are 

2 /-^/^ 

-L/2 

From Eqs. ([5]) and ([H]) we see that these coefficients are the solution to the initial-value 
problem 

^ = --(y - -')^ + If^'y^ - \ [(2- + ^)-ILfySn. (9) 

with initial condition 

^„(l/,0)=7n, (10) 

where 

Defining a new time scale r and a new volatility variable v by the change of scale 

r = at, v = {2a/P)y, (12) 

the problem above reads 

§^ = -(. - ^)f^ + .^ - (/3./2L)^.5., (13) 

Sn{v,0) = jn, (14) 

/3„ = (A;/a)(2n + l)7r, (15) 

^ = {2a/e)m^. (16) 



and 



where 



and 



Before proceeding further let us remark that the parameter 6, which turns out to be 
crucial for the escape problem at hand, can be regarded as the "dimensionless normal level" 
of volatility. It represents a balance between the tendency toward the normal level measured 
by am^ and the volatility fluctuations quantifled by k"^ (see discussion in Sec. lIVBI regarding 
the cases 9 < 1 and ^ > 1). 



The problem posed by Eqs. (TT^ - (IT^ is solved by the function 

Sn{v,T) =7„exp{-A„(r) - B„(r)t;} , (17) 

where A„(r) and Bn{T) are functions of time to be determined. In effect, plugging it into 
Eq. flT^ we see that Eq. flT7|) is the solution to the problem provided that 



A„(r) = 9 / B^{s)ds, (18) 

and Bnir) obeys the Riccatti equation 

5„ = -5„ - El + (/3„/2L)^ (19) 

with initial condition -B„(0) = 0. 
In the Appendix |X] we show that 



and 



jji^T + In 



An 



(20) 



where 

A„ = ^/iTWlW. /^± = (An ± l)/2. (22) 

Therefore the solution to the escape problem for the two-dimensional Heston SV model 
is 

oo 

S{x, f , r) = ^ 7„ exp {-A„(r) - i?„(r)t;} cos [{2n + 1) vrx/L] . (23) 

n=0 

Figure [T] shows in a three-dimensional plot this SP as a function of the return x and volatility 
variable v for r = 0.1 |l9l |. 

In the asymptotic regime, either for long or short times, the SP is somewhat simpler. 
Thus when r ^ 1 (i.e., t ^ a^^) we have Anij) ~ On^r and Bnij) ~ /i_. Hence 

oo 

5(a;,i;,r) ~^e-'^-(^"+^)cos[(2n + l)7ra;/L], (r > 1). (24) 

n=0 

On the other hand for short times r -C 1 (i.e., t <^ a^^) we write e^'^"* = 1 — A„t + O(t^) 
and taking into account that /i_ + /i+ = A„ and /i-/i+ = — (/3n/2-L)^, we see from Eqs ( 120|1 
and (EU) that 

A„ = e[^^^r + ln(l - /._r)] + 0{r^), i?„ = i^^iZ^^ + 0(r2), 

i — /i_r 
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FIG. 1: (Color online) The survival probability S{x,v,t) given by Eq. ()23p with r = 0.1 (t = 
2.22 days) and L = 0.01 in terms of return x and volatility v. Parameters of the model: 9 = 1.25, 
a = 0.045 day"\ m = 0.093 day"^/^ and k = 0.0014 day'^ Recah Eq. ([16]) and note that there 
only exist three independent parameters. 



whence 

S{x,v,t) 



OO -. 

n=0 ^ ^ 



^_T) 



{l3n/2Lfv ' 
1 — fj,_r 



exp<;- g/i, + ^7"^ ".^^_" T ^ cos [{2n+ I) Tcx/L], (r < 1). 

(25) 
In Fig. [2] we represent the exact SP, Eq. fl23|) . at x = and for fixed times, in terms 
of the volatility. The plots confirm , as hinted by Eqs. (1241) and (1251) . that the SP decays 
exponentially with the volatility for both short and long times. The characteristic exponent 
of this decay depends on the value of 6, being larger for smaller 6, i.e., larger k (cf. Eq. 
(IT6l)). Moreover, as r ^ 1, when the volatility is small, the higher survival probability 
corresponds to the case when 6 is smaller. This is a distinct behavior with respect to the 
rest of situations. 




FIG. 2: (Color online) The survival probability S{x,v,t) given by Eq. (I23p at x = as a function 
of the volatility v. Left plot shows the case when r = 0.1 (t = 2.22 days) and L = 0.01. The 
figure on the right exhibits the case when r = 100 and L = 0.1. The straight lines correspond 
to the exponential decay with v mentioned in the main text. Parameters of the model are a = 
0.045 day" , m = 0.093 day" ' and the three different values of the parameter 6 = {2a/k'^)m? 
provide three different values for k accordingly. 
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FIG. 3: (Color online) The mean-escape time T{x,v) given by Eq. (I27p in terms of return x and 
volatility variable v. Parameters of the model: 6 = 1.25, a = 0.045 day" , m = 0.093 day" ' and 
A; = 0.0014 day"^ 



III. THE MEAN ESCAPE TIME 



The survival probability S{x, v, t) provides maximal information on the escape problem of 
the two-dimensional process (X(t),F(t)). Indeed, the probability density function f{t\x,v) 



of the escape time is related to the SP by [iSj 



dS{x,v,t) 
f{t\x,v) = , 

and all moments of the escape time can be obtained through the SP. Thus, for instance, the 
mean escape (or exit) time (MET) is given by 

T{x,v) = / S{x,v,t)dt. (26) 



Jo 
For the Heston model we see from Eq. ([23]) that the two-dimensional MET, T{x,v), can 
written in terms of the following Fourier series 



^ oo 

T{x, t;) = - y T^{v) cos [{2n + 1) ttx/L] , (27) 

where 



a 

n=0 



Tn{v) = 7n / exp{-v4„(r) - Bn{T)v}dT. 
Jo 

Using Eqs. ( l20l) - (!2Tl) and some simple manipulations, which involve the change of variable 



e '^"^, yield 








Tn{v) = 


" ^^% Jo [1 + (/^-^+)er """n 


-/i+ 


[ i-e 1 


_i + (/i_//i+)e. 



(28) 

Figure [3] provides a three-dimensional representation of T{x,v) based on the numerical 
computation of Eqs. ( I27j) - (!28l) [20]. A noticeable aspect worth stressing is shown in Fig. H] 
where we two projections of the MET T{x, v) are depicted either for small and large volatility 
and also for three different values of the normal level 6. Thus when volatility is very large, 
V = 1300, the left plot in Fig. H] shows that the larger 6 corresponds to the longer MET. In 
the opposite case of very low volatility, v = 0.001, the right plot shows that this behavior is 
reversed, for now 9 = 1 corresponds to a longer T{x,v). This anomaly is also observed in 
Fig. [5] when v < 0.1. 

Having obtained the expression for T(x, v) as given by Eqs. ( 1271) and (l28l) . let us proceed 
to elucidate how is the dependence of the MET on the volatility. This is a meaningful 
question from a practical point of view, for market behavior depends critically on volatility. 
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FIG. 4: (Color online) The mean-escape time T(x, v) given by Eq. (I27p as a function of the starting 
return x. Left plot shows the case when v is large showing a perfect hierarchy where larger d means 
larger MET. Right plot shows how the 9 = 1 case breaks this hierarchical order for small enough 
values of v. Parameters of the model are a = 0.045 day" , m = 0.093 day" ' and the three 
different values of the parameter 9 = {2a/k'^)m'^ provides three different values for k accordingly. 
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FIG. 5: (Color online) The mean-escape time T{x,v) given by Eq. ()27p at x = as a function 
of the starting volatility variable v. The drawings illustrate that the MET saturates at certain 
maximum value when v = 0. On the other hand, the straight line clearly shows that the MET 
decays as 1/v with increasing volatility as we prove in the main text. Parameters of the model are 
the same than those of Fig. HI 
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Intuition tells us that the escape time must tend to zero as the volatility increases and a 
quick glance both at Eq. fl^Hj) and Fig. [S] confirms this, but how is the form of this decrease? 
On the other hand the behavior of the escape time if the volatility is low is also relevant, will 
T{x, v) grow without bound as t; — i> 0? or will it tend to a finite, albeit maximum, value? 
We will next answer these questions. 

Let us first obtain the behavior of the MET when f = 0. In such a case Eq. fl28p reads 

A^-i n t-i+M-e/A„ 



and using the following integral representation of the Gauss hypergeometric function 2l| : 

i^(«, h- c; z) = ^^^^^^^l_ ^^ f^ e-\l - 0'-'-\l - ^zr'^di {c>b>0), (29) 

we have 

We therefore see that the mean escape time T{x, v) tends to a finite quantity when v = 0. 
Let us now turn to the case of increasing volatility. In this situation it is convenient to 
perform the following change of integration variable in Eq. f l2H|) : 

then (recall that /i+/i_ = (/3n/2L)^) 

Jo 



where 



g{z) = (1 - /i+^)-^+'^-^/^"(l + /i_^)-i+^-^'-^/^^ 



As f ^ oo the exponential term falls off quickly and we may safely change the upper 



integration limit l/yU+ by cxd. Using then Watson's lemma we write 22 | 



oo 



(2L//3, 



\2fc 



r..(.)~%E9'"(0)^^^P^- (31) 

k=o 

Up to the leading order ((7(0) = 1) 

Tn{v) ~ 7n {2Llf3nf {l/v) + O {l / V^) , 
11 



or (cf. Eqs. (|TT|) and (fTUj)) 



Tn{v 

Therefore 



~^5^(^rTiF<^''"'^°(^''"^)- ('^' 



n=0 ^ '^ 



The series on the right can be summed with the result 23 1 

nx,v)r^—[{L/2f-x']+0{l/v'). (33) 

This is a remarkable result since shows that for large volatility the MET has the same 
form than that of the Wiener process (see Sect. lIVp . Moreover the two-dimensional MET 
decreases linearly as 1/v. 

The behavior of T(x, f) with volatility is clearly seen from the numerical evaluation of 
the exact MET given by Eqs. flTTI) and (1^ . Figure shows, in a log-log scale, how the 
MET saturates to a maximum when v tends to zero while for large volatility T(x, v) is well 
fitted with a power law with exponent -1 which confirms the asymptotic expression (l33l) . 

IV. AVERAGING THE VOLATILITY 

In real financial data the volatility is, in fact, a hidden variable which has to be measured 
in an indirect way 2S[ ■ It is therefore of great significance to know whether the price of an 
asset remains inside a given interval regardless its volatility. In physics the analog to this 
question would be knowing the survival probability for the position of a Brownian particle 
without worrying about its velocity [7]. Considering the entanglement between return and 
volatility (or position and velocity) this is certainly a difficult question and one often has to 
rely on approximate answers. Fortunately the latter is not the case in the Heston model as 
we shall see next. 

A. The survival probability of the return 

In order to obtain the SP of the return, S{x, r), regardless the value of the volatility, we 
have to average the volatility away from S{x, v, r). We will do this by assuming that, at the 

12 
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FIG. 6: (Color online) The survival probability of the return, S{x,t), given by Eq. (I4ip with 
L = 0.01 as a function of the return x and time r = at. Parameters of the model: 9 = 1.25, 
a = 0.045 day-\ m = 0.093 day'^/^ and k = 0.0014 day'^ 



time we measure the return, the volatihty process has reached the stationary state 25|. We 
therefore define S{x, r) as the average: 



S{x,t) 



S{x,v,T)pst{v)dv, 



(34) 



where Pstiv) is the stationary probabihty density of the volatility. For the Hestson model 
this density is the normalized solution of the Fokker-Planck equation 



d 
dv 



, /IN ^ 

[v-B)^ —V 
dv 



Pst{v) = 0, 



which is given by the Gamma distribution: 



Pst[V) 



no) 



v'-^e-\ 



(35) 



Note that Q is the stationary variance of the volatility variable v. Indeed, from Eq. fl35|) 
we see at once that Q = ('y^)st — ('^)sf Observe also the changing shape of the stationary 
distribution (specially as f ^ 0) according to whether 6' < 1 or 6* > 1; a fact that, as we 
shall see below, has consequences on the behavior of the MET. 
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From Eqs. ([7]) and (IMI) we get 

oo 

S{x, t) = J2 Sn{r) COS [{2n + 1) nx/L] , (36) 

n=0 

where 

Sn{r) = / Sn{v,T)pst{v)dv, 

Jo 
which, after making use of Eqs. ( fTTl) and ( l35l) . yields 

We will write this Fourier coefficient in a more convenient form. Let us first note that by 
applying Eq. (I20l) we can write 

e-^"(-)=f ^""'"'^^ y. (38) 



On the other hand, from Eq. ( l2Ti) we see 

/i+(l + /i-)+/i-(l-;^+)e-^"" 
/i+ + /i_e ^"^ 

but 1 + /i_ = /i+ and 1 — yU+ = — /i_ (cf. Eq. f l22|) ). Hence 



1 + B„ r = ^+ /" ■ . 39 



Plugging Eqs. ([321)-® into Eq. §7i) we have 



A„e-'^- 



^n(^)=7..(^ ^^_^,_^_^„J ' (4°) 



and therefore 



°° /A -/t-T \ ^ 

^(x, r) = 5^ 7n . _"' „^^, cos [(2n + 1) nx/L] , (41) 

„=o V^+ ^-^ / 

which constitutes the exact expression for the SP of the return. 

We will now show the asymptotic time behavior of S{x, r). We easily see from Eq. (HTl) 

that for long times, r ^ 1, the asymptotic form of the SP is 

S{x,T)-J2^n(^) e-^'^-"cos[(2n + l)7rx/L], (r > 1), 

n=o ^^+/ 

while for short times r -C 1 and after taking into account (cf. Eq. (!22l) ) 

/i^ - /i'_e-^'- = A„(l + /iV) + 0(r2), 
14 



(42) 
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FIG. 7: (Color online) The survival probability S{x,t) given by Eq. (I4ip with L = 0.01 as a 
function of the time t = at and starting from the return midpoint x = 0. Parameters of the model 
are: a = 0.045 day" and m = 0.093 day" ' . Notice that since a = 0.045 day" then r = 1 
corresponds to an actual time of t ~ 22 days. The straight lines plotted on the lowest curve are 
the asymptotic approximations given by Eqs. (Il2|l (r 3> 1) and (fi3ll (r ^ 1). 



we get 



S{x,t) -X]^"7f 



-Bll-T 



COS [{2n + 1) TTx/L] 



(r « i: 



n=0 



The approximate expressions for S{x, r) given in Eqs. (I42l) - fj43l) suggest an exponential 
decay (essentially governed by the normal level 6) either for short and long times. This is 
confirmed by the numerical evaluation of the exact SP, Eq. (1411) . which we present in Fig. 
[71 We clearly see there two different exponential decays which match those shown in Eqs. 
dH-dH. 



B. The mean escape time of the return 



In terms of the survival probability S{x, t) the mean escape time is given by 



T{x) = I S{x,t)dt. 
'o 
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Combining this equation with Eq. (HTl) we see that T{x) is written as a Fourier series of the 
form 

^ oo 

T{x) = -y^TnCos[{2n + l)nx/L], (44) 



a 

n=0 



with Fourier coefficients given by 



The integral appearing in the right hand side of this equation is evaluated by performing 
the change of variables ^ = e"^""^. We have 

l^+ Jo 
and taking into account the integral representation of the Gauss hypergeometric function 
given in Eq. ( l29l) we get 

Finally, the MET is given by (see Fig. [H]) 

n^)-ig^(|)"K''^^^+^^S)~^|''"+^''^/'|' <'"' 

From a practical point of view an interesting property to look at is the behavior of the 
MET as a function of the span L specially for short and large values of L, being the latter 
closely related to financial defaults or uprisings depending on the sign of x. We will thus 
consider the two limiting cases: (a) L ^ and (b) L ^ oo. 

(a) In the case of small span the Taylor expansion as L — >■ of Eq. (HUj) leads to the 
following asymptotic expression (see Appendix [B] for details) 



T(x) 



\e+i 


9<1, 


-L'^lnL 


9 = 1, 


L' 


e>i. 



{L -. 0) (47) 



We see that in this case the behavior of the MET is governed by the (dimensionless) normal 
level 6 which coincides with the stationary variance of the volatility variable v. Let us recall 
that a similar situation arises for the stationary distribution since, as seen in Eq. 
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FIG. 8: (Color online) The mean-escape time T{x) based on the exact expression ()46p as a function 
of the starting return x. The left figure shows the case when L = 0.1 while the right plot shows 
the case when L = 0.01. In both cases the larger the the lower the MET. We also draw the MET 
corresponding to the Wiener process (note that the latter is always shorter than Heston's MET). 
Parameters of the model are the same than those of Fig. [H 

Pst('y) behaves in a different way according to whether the normal level is greater or lower 
than 1. Note that (cf. Eq. (TTUj) ) Q < \ implies m^ < k'^/a, that is, volatility fluctuations 
-represented by the vol-of-vol k- are wilder than the tendency toward the normal level given 
by am'^. On the other hand, when this tendency is greater than the volatility fluctuations 
{6 > 1) the MET grows quadratically with L independent of the normal level m but with 
slope depending on the vol-of-vol through the combination k'^/a. All of this is exemplified 
in Figs. [UKTUl where we plot, based on the exact expression (HUI) . the MET as a function of 
the span L. 

(b) Let us now look at the behavior of the MET with increasing span. Unfortunately 
this case is more difficult to deal with since L appears in the Fourier series solution basically 
through the combination {2n+l)/L and any effect due to L ^ oo is neutralized by increasing 
values of n which, in turn, are needed to sum the Fourier series. The case (a) above turns 
out to be workable because the limits L ^ and n ^ oo are compatible. 

On the other hand, numerical calculations of the exact MET fHUl) shown in Figs. [UlITU] 
clearly indicate that the MET grows quadratically with the span regardless the value of the 
normal level m: 

T{x) ~ L\ {L -^ oo). (48) 

We remind that we have already encountered this behavior at the end of Sect. IIIII when 
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analyzing the two-dimensional MET, T{x,v), for large volatility (cf. Eq. (1331) ). In the 
Appendix O we justify this quadratic growth by means of a heuristic argument. 

We also note that now, contrary to the case of small span, the slope is independent of 
k'^/a and all the cases which have the same m merge into a single curve (cf. Figs. iQllTOl) 

C. The Wiener process 

For many years the most ubiquitous market model has been the geometric Brownian 



motion which was proposed by Osborne in 1959 26|. In this model the price P{t) obeys the 

stochastic differential equation 

dP(t) 

—i-i = udt + adW{t), 

r[t) 

where z/ is a constant drift, a is the volatility (a constant as well) and W{t) is the Wiener 
process. In terms of the zero- mean return X{t) defined in Eq. ([T]), the model reads 

dX{t) = adW{t). 

In other words, X{t) is the Wiener process with variance cr^. 

In view of the widespread use of this market model among practitioners and even aca- 



demicians 271] , we find it convenient to compare the findings on the escape problem of the 
return discussed in this section with those of the Wiener process. This, in turn, may pro- 
vide a test on the appropriateness of the assumption of stochastic volatility for real market 
models. 

Let us thus suppose that the zero-mean return is described by the Wiener process and 
denote by So{x,t) its survival probability inside the interval — L/2 < X(t) < L/2. This 



function obeys the equation [18 1 



a5o _ 1 2^ 

with initial and boundary conditions 

5o(x,0) = l, Soi±L/2,t)=0. 

Proceeding as we have done before we look for a solution to this problem in terms of a 
Fourier series. In this way one easily obtains 

So{x,t) = -y^i, ^exp{-[7rL(T(2n+l)]V2}cos[(2n+l)7rx/L], (49) 

71 ^ — ^ 2n + 1 

n=0 
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FIG. 11: (Color online) Return survival probability as a function of the scaled starting return x/L 
when L = 0.05. We represent the S{x,T/a = 5 days) of the Heston model given by Eq. (I34p for 
several values of 6. We also plot the SP corresponding to the Wiener process, 50(0, t = 5 days), 
given in Eq. (j49|) . Parameters of the Heston model are a = 0.045 day~ and m = 0.093 day^ ' . 
For the Wiener case we suppose that volatility is equal to the normal level a = m. 



and the MET is 



Ux) = -[{L/2r - x' 



o^ 



(50) 



In Fig. [TT] we plot the Sq{x^ t) given by Eq. (H9l) in terms of the return x and for a fixed 
time t = 5 days. In Fig. [12] we do the same but as a function of time and for a fixed return 
X = 0. In the both figures we also represent the Heston's SP, S{x,t), given in Eq. fj4T]) . We 
see that the survival probability is always higher under stochastic volatility than when the 
volatility is constant; although for a greater normal level 6 this difference becomes smaller. 

Thus, for instance (cf. Fig. [T2l) when 6 = 1.25 the survival probability of the Wiener 
process in one day starting at x = with span L = 0.01 is just 6*0(0, 1 day) = 0.026; for 
the Heston model this probability is 8 times higher: 5(0, 1 day) = 0.208. In two days the 
difference is even higher: S'o(0,2 days) = 0.0005 versus ^(O, 2 days) = 0.095. 

This difference is also detected in the MET. Thus in Fig. [S] we see that the Heston's 
MET is invariably longer than that of Wiener's. In other words, the Wiener process escapes 
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FIG. 12: (Color online) Return survival probability as a function of the scaled time r. We represent 
the decay 5(0, at) with time for the Heston model given in Eq. ()34p when 9 = 1.25 in comparison 
with the Wiener model SP, 5*0(0, t), provided by Eq. (j49p . In both cases we assume that starting 
return is the midpoint of the interval L = 0.01. Parameters of the model are the same than those 
of Fig. H 

faster than the Heston SV model. 

Therefore, and contrary to intuition, the assumption of stochastic volatility, notwith- 
standing occasional bursts, seems to stabilize prices after a certain number of time steps. 

V. SUMMARY AND CONCLUSIONS 

We have studied the escape problem of the return under the assumption of a stochastic 
volatility given by the Heston model. The problem is fully characterized by knowing the 
survival probability, S{x,y,t), of the bidimensional process {X{t),Y{t)) inside the strip: 

-L/2 < X{t) < L/2, < Y{t) < oo, 

where X{t) is the zero-mean return and Y{t) is the volatility variable which, for the Heston 
model, is related to the volatility a by a(t) = ^^/Y{t). The survival probability obeys the 



21 



backward Fokker-Planck equation (^ with initial and boundary conditions given in Eq. ©• 
We have been able to exactly solve the problem by means of the Fourier series expansion 
given in Eq. (l23l) . 

Once we have the solution for S{x, y, t) another interesting and most useful quantity to 
know is the mean escape time. For the entire process given by the return and the volatility, 
the MET T{x,y) is also exactly given by a Fourier series (cf. Eqs. fl27l) - (!2H|) ). We have 
shown that, as the volatility decreases, T{x, y) tends toward a maximum, albeit finite, value. 
Moreover, as the volatility increases the MET decreases following the hyperbola 1/y which 
is quite remarkable because this is exactly the behavior of the MET with the volatility had 
the retun followed the Wiener process (i.e., constant volatility) instead of the Heston model. 

Real financial data consist of time series of prices and volatility is not directly recorded 
and only observed in an indirect way. This hidden character makes it worth averaging out 
volatihty from the expressions of S'(x, y, t) and T(x, y) and thus solving the escape problem 
for the return alone. The assumption to be made is that the volatility has reached the 
stationary state; in the Heston model the latter is characterized by the Gamma distribution, 
Eq. (ESD. 

Following this way we have obtained exact expressions of S{x, t), Eq. fj4T|) . and T(x), Eq. 
(HUj) . both in terms of Fourier series. The SP has two different exponential decays: one for 
long times and another, which is faster, for short times. We have been able to get analytical 
expressions of both decays. 

We have next analyzed the behavior of T{x) as a function of the span L, specially for 
short and large values of L. The latter case is specially significant because large values of 
L are associated to financial uprisings or defaults. We have shown that the behavior of the 
MET as L ^ depends on the normal level 9 and it is given in Eq. fHTj) . On the other 
hand, when L ^ oo the MET grows as L^ independently of the normal level. 

Therefore, when ^ < 1 (i.e., if volatility fluctuations are greater than the tendency toward 
the normal level) we have a "crossover" in the MET, from L ^ to L —>■ oo, of the form: 

T{x) ~ L'+^ -^ T{x) ~ L^, {6 < 1). 

On the other hand for > 1 (the tendency to relax toward the normal level is now stronger 
than the fluctuations of the volatility) there is no such crossover, since T{x) ~ L"^ for both 
small and large values of the span. Again, this quadratic dependence is the same as if the 
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return would had been described by the ordinary Wiener process. 

We have finally compared the return SP and MET to those provided by the Wiener 
process. In other words, we have confronted the escape problem of the Heston model with 
that of the Wiener process. Our main finding is that Heston's SP is bigger and Heston's 
MET is longer than those corresponding to the Wiener process. This, at first sight, is 
contrary to intuition because a random volatility, despite occasional bursts, would seem to 
stabilize prices to a greater extend than a constant volatility. However, let us recall that in 
S{x, t) and T{x) the volatility has been averaged around its mean value that is precisely the 
normal level 6 and, if 6 is not very large, this fact may be the responsible for the stabilization 
of the return. 

A final remark. The Heston model is one among several possible candidates aimed to 
describe a realistic price dynamics. The question of which SV model is more appropriate as 
a market model is still an open question J28|. We have chosen the Heston model to carry 
on the present development because, as we have seen, it allows for an exact treatment. In 
a forthcoming work we will present an approximation scheme in order to study the escape 
problem for a wider class of models. 
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APPENDIX A: FUNCTIONS ^„,(r) AND 5„(r) 

To obtain the functions A„(r) and Bn^r) we must solve the Riccatti equation 

B„ = -Bn -Bl + (/3n/2L)^ (Al) 

with initial condition -8^(0) = 0. To this end we define a new function Znir) related to 
Bn{r) by 

Bn{r) = |. 
then Z{t) obeys the linear equation 

Z + Z- {(3j2LfZ = 0, 
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whose solution reads 

where Ci and C2 are arbitrary constants and 



/i± = (A, ± l)/2, A„ = Vl + 4/32 
The expression for Bnir) is thus given by 

and the initial condition Bn{0) = yields 

C2/C1 = Ai-//i+. 



Hence, 



Bn{r) = /i_ 



1-e 



-A„r 



(A2) 



l + (/x_//.+)e-^"-' 
which is Eq. f ETjl . 

Plugging Eq. ( lA2p into Eq. (fTSl) and setting ^ = e~^"'^ as a new integration variable we 

get 

1-e 



^n(r) = ^ /' TT^-r-'r-^d^ 



but 



1-e 



Hence (recall that fi^^ + fi^ = A„) 



A„ y,-A„. e[i + (/i-//^+)e] 

rfe = Ine - (1 + /^+//i-) ln[l + {fi-/f^+n 



yU_r + In 



aT 



'A„T 



which is Eq. (EO 



(A3) 



APPENDIX B: BEHAVIOR OF THE MET FOR SMALL SPANS 



Let us suppose L — >• 0. From Eq. fl22l) we see that 



A„ = i/3jL) [I + {L/l3nf/2 + 0(L^)] 



and 



/i± = {(3J2L) [I ± (L//3„) + 0{L' 
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Hence 

,2 



2 -1-(4L//3„) + 0(L2), (B1) 



and 

1 ( ^y 
72 



/i 

^ = [1 - (L//50 + 0(L2)]/2, (B2) 



22^+3(L//3„)^+i [1 + (1 - 2e){L/Pn) + 0(L2)] . (B3) 



(i) Suppose that 6 <1. Using Eq. (JBip we write 






fl^.^ii + ^i^l^fl^.^ii + ^nl+oiD, 



but [211 



F(a,6;c;l) = ^^|^^^\-^, (c-6-a>0). (B4) 

7(c — a)T[c — b) 

Note that condition c - 6 - a > imphes 61 < 1. We thus find (cf. Eqs. flBT]) -flB2l)) 

^ ^^' A^' ' + AT' ^4 J = r(i-g/2) '^ + ^(^)]' (^ < ')• (^') 

Plugging Eqs. ( ]B3P - (1B5I) into Eq. ( l46l) and taking into account Eqs. (ITTl) and (iT6l) we finally 
get 

n=0 

where 



1—1)" 



/2^^+^\ /2ay+^ r(i + g/2)r(i-g) 
^^ = [-^^J UJ r(i-0/2) — ■ ^^^^ 

For a; = we have 

T(0)=iriL''+^[l + O(L)], (^<1), (B8) 

where 

oo 
n=0 

and for small values of the span L the MET grows as a power law with exponent related to 
the normal level of the volatility 9 < 1. 

(ii) Suppose now that 6 > 1. We employ the following property of the hypergeometric 



function 21 1 



F{a, b; c; z) = {1 - zy-''-'F{c -a,c-b;c; z) 
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and write 






1- 






2 \ 1-fi 



F(l-^ + ^a;^^--^- 



' ' A ' ,,2 



which, after using Eq. fIBip . reads 



F 1 



A„ yu^ 



^■^^t^^l-°<-) 



Note that we can apply Eq. ( JB4I) since condition c — a — 6 > now iniphes ^ > 1 . Hence 
(cf. Eq. (12}) 









1-e r 



9-1 



r(l + 0/2) + O(L) 



Substituting this into Eq. (H6ll and taking into account Eq. (1B3I) we obtain 



T(x)=iV2L2^ 



n=0 



(2n + 1) 



cos[(2n + l)7rx/L][l + 0(L)], (^ > 1) 



where 



No 



2' \ f2ayT{l + 9/2) 



For a; = we have 



where 



^naO J \nkj 1 — 9 
T{0) = K2L^[l + O{L)], {9>1), 

oo 
n=0 



(B9) 



(BIO) 



(Bll) 



Therefore, in this case the average escape time grows quadratically with the span -as if the 
zero-mean return would have followed the simple Brownian motion- independently of the 
value of the normal level of volatility 9 > 1. 

iii) When ^ = 1 we utilize the following series expansion of the hypergeometric function 



m 



F(a, b:a + b: z) 



^^''^^^E%^]2v^(^+i)-v^(«+^)-^(^+^)-i^|i-^i](i-^r' 



r(a)r(6)— m 



which, when z —>■ 1, yields the following approximation 

T{a + b) 



F(a, b;a + b; z) 



V{a)V{b) 



lull -z\ + 0(1). 
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Hence, as L — i> 0, 

,2 



F(l,^;l + ^;^)=-(l/2)ln(4L//3„) + 0(LlnL), 



whence 



T{x) = ^—^ y2^[- ln(4L//?„) + 0(L In L)] cos [(2n + 1) vrx/L] , (^ = 1), (B12) 

and there is a logarithmic growth with the span when 6 = 1. 

The resuhs above (cf. Eqs. ( ]B6p . ( 1B9I) and ( 1B12I) ) are summarized in Eq. ( 1471) . 



APPENDIX C: BEHAVIOR OF THE MET FOR LARGE SPANS 

Let T(x, y) be the MET of the joint process (X(t), F(t)) out of the strip -L/2 < X{t) < 
L/2, < Y{t) < oo. In terms of SP S{x,y,t) the MET is given by 

/■oo 

Tix,y) = / S{x,y,t)dt. 
Jo 

From Eqs. (I5])-(J6]) we easily see that T{x,y) is the solution to the boundary-value problem 

T(±L/2,|/,t)=0. (C2) 

Now our heuristic argument: large values of the span L are equivalent to small values of x, 
but in this situation (as long as y is not too small) d'^T/dx'^ is the dominant term in the left 
hand side of Eq. ( JCll) . This allows us to approximate the MET T{x,y) ~ To{x,y), where 
the "outer" approximation [29] To{x,y) is the solution to 

ly^ = -^^ n±L/2,y,t) = 0. 



That is, 



and the MET grows as 



To{x,y) = -[{L/2y-x'], 



T{x,y)r^L'^ (L^oo). (C3) 
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